Charm meson resonances and D — > V semileptonic form factors 
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Using limits of large energy effective theory and heavy quark effective theory we propose a simple 
parametrization of the heavy to light H — -> V semileptonic form factors. Then we reconsider D — > 
Vlvi decays within a model which combines heavy meson and chiral symmetry. In our Lagrangians 
we include contributions coming from excited charm meson states, some of them recently observed. 
Within this framework we determine all parameters describing the shapes of the form factors and 
calculate branching ratios and helicity ratios for all D — > Vive, decays. 
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I. INTRODUCTION 

Presently, one of the most important issues in hadronic 
physics is the extraction of the CKM parameters from 
exclusive decays. An essential ingredient in this approach 
is the knowledge of the form factors' shapes in heavy to 
light weak transitions. Usually, the attention has been 
devoted to B decays and the determination of the phase 
of the V u b CKM matrix element. At the same time in 
the charm sector, the most accurate determination of the 
size of V cs and V c d matrix elements is not from a direct 
measurement, mainly due to theoretical uncertainties in 
the calculations of the relevant form factors' shapes. 

Recently CLEO and FOCUS have published interest- 
ing results on D° — > it~(. + vi and D° — ► K~t r Vt de- 
cays pi Q Their studies indicate that the single pole 
parametrization of the relevant form factor cannot ex- 
plain their data very well, leading to unphysical pole 
masses. Both experimental groups have also attempted a 
modified pole fit, which was first put forward for B decays 
in Q , and their results suggest the existence of contribu- 
tions beyond lowest lying charm meson resonances. 

On the other hand we have recently |4| reconsidered 
D — > Ptvi decay form factors within a framework which 
combines heavy meson and chiral symmetries (HM;yT) 
and includes in the interacting Lagrangian contributions 
coming from excited charm meson states. We have found 
that a two-poles shape of the relevant form factor can 
be successfully accommodated within HM^T when ex- 
cited meson states are included into the model. In our 
approach the first pole is described by the lowest lying 
vector resonance, as in the original idea , while for the 
second one we assume complete saturation by the next 
vector state which we include in our Lagrangian. In do- 
ing this we anticipate the discrepancies from the general 
two-poles procedure, in which the second effective pole 
should account for all other excitations that might be ex- 
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changed in the t-channel, to be small and encoded in the 
parameters of the model. The unknown parameters have 
been obtained by fitting the experimental results for the 
branching ratio. The assumed pole behavior agrees well 
with experimental results confirming our anticipation of 
small saturation error. 

In addition to studies of heavy to light pseudoscalar 
meson weak transitions (H — ► P), transitions of heavy 
pseudoscalar mesons to light vector mesons (H — > V) 
such as D s — > (fiive an D s — > K*lvg offer an opportunity 
to extract the size of the relevant CKM matrix elements. 
We continue with our study and re-investigate vector and 
axial-vector form factors in D — » Vivi decays within a 
similar framework as in the case of D — > P£ff. The 
H — ► V transitions were already carefully investigated 
within many different frameworks such as perturbative 
QCD PJI, QCD sum rules 0, B H E El E3 , lattice 
QCD Ellll 13 111 , a few attempts to use combined 
heavy meson and chiral La gran gians (HMvT) 0, Il8j . 
quark models fljj 12(1 I2TL l22j. large energy effective 
theory (LEET) |23l and soft collinear effective theory 
(SCET) HEllMElllllli. Each of these approaches 
has only a limited range of validity. For example, the 
QCD sum rules, LEET and SCET are suitable only for 
the low q 2 region while lattice QCD and HM;yT are suc- 
cessful for maximal q 2 . It is important to note, that 
currently lattice QCD and QCD sum rules are the only 
approaches that enable the computation of form factors 
solely from first principles. On the other hand quark 
models usually involve parameters which have little phys- 
ical correspondence to the underlying theory of QCD. 

The experimental situation in D — > Vlvi has not 
changed a lot in the last few years, but recently it has 
been gaining pace yfl HU, H2, H3, H3, |35| , and hopefully 
more results on the q 2 shape of the form factors will be 
available soon. Unlike in the case of H — > P weak tran- 
sitions, no general parametrization of the form factors, 
relevant to H ~ * V weak decays has yet been proposed. 
Usually a simple pole behavior of all the form factors is 
assumed when extracting values of the form factors at 
q 2 = from experiment or extrapolating results of dif- 
ferent theoretical approaches. 
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Recently, the spectroscopy of charm mesons has been 
enriched by discoveries of many new charm meson res- 
onances. BaBar 36j collaboration has announced a 
new, narrow meson D s j(2317) + . This was confirmed 
by Focus 33 and CLEO [38| which also noticed an- 
other narrow state, £> s ,/(2463) + . Both states were con- 
firmed by Belle [3!j. Finally, Selex |4jj nas announced 
a new, surprisingly narrow state -D^j (2632). The states 
D s j(2317) + and D s j(2463) + are already being identi- 
fied to belong to the (0 + , l + )spin-parity doublet of the 
D s mesons while the D^j (2632) state has been proposed 
as the first radial excitation of the Z?*(2112) with the 
spin parity assignment 1~ [ill E2I |43| . 

The purpose of this study is to (1) devise a general 
parametrization of all the form factors relevant to H — > V 
weak transitions which would take into account known 
experimental results on heavy meson resonances as well 
as known theoretical limits of heavy quark effective the- 
ory (HQET) and LEET relevant to H -> V weak tran- 
sitions; to (2) investigate contributions of the newly dis- 
covered charm mesons to D — > V semileptonic decays 
within an effective model based on HM^T by incorpo- 
rating the newly discovered heavy meson fields into the 
HM^T Lagrangian and utilizing the general form factor 
parametrization. We restrain our discussion to the lead- 
ing chiral and l/mn terms in the expansion, but we hope 
to capture the main physical features about the impact 
of the nearest poles in the i-channel to the q 2 -dependencc 
of the form factors. 

In Sec. n we revise the common weak current matrix 
element decomposition relevant to transitions between 



pseudoscalar and vector mesons and introduce a form fac- 
tor decomposition, which is independent of the mass of 
the pseudoscalar meson and thus convenient for studying 
H — > V weak transitions. In Sec. Ill we derive a general 
H — ► V form factor parametrization drawing from both 
known experimental properties of heavy mesons as well 
as from known theoretical scaling laws and form factor 
relations in the limit of the infinite heavy meson mass. 
Sec. IV describes the framework we use in our HM^T 
calculations: we write down the HM^T Lagrangian for 
for heavy and light mesons and extend it to incorporate 
new heavy meson fields. In Sec. V we calculate the val- 
ues of the D — > V semileptonic form factors near zero 
recoil within HM^T and extrapolate our results to larger 
recoils using the general parametrization of Sec. Ill and 
by saturating the effective poles with physical masses of 
experimentally known or theoretically predicted charmed 
resonances. Finally, a short summary of the results and 
comparison with other approaches as well as with exist- 
ing experimental data is given in Sec. VI. 



II. PARAMETRIZATION OF H -» V CURRENT 
MATRIX ELEMENT 

A frequently encountered decomposition of the current 
matrix elements relevant to semileptonic decays between 
a heavy pseudoscalar meson state \H{pu)} with momen- 
tum p v H and a light vector meson state | V(pv, ey)) with 
momentum p v and polarization vector e v is 



(V(e v ,p v )\q^Q\H(p H )) 
(V(e v , Pv )\q^ 5 Q\H(p H )) 



2V(q 2 



ma + my 
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ty ■ q 
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Al(q 2 ) 



(1) 



where q — u, d or s are the light quark fields, Q = b or c 
denote the heavy quark fields, q v = {pn — PvY is the 
exchanged momentum and q 2 is the exchanged momen- 
tum squared. Here V denotes the vector form factor and 
is expected to be dominated by vector meson resonance 
exchange, the axial A\ and A2 form factors are expected 
to be dominated by axial resonances, while Aq denotes 
the pseudoscalar form factor and is expected to be dom- 
inated by pseudoscalar meson resonance exchange |l9j |. 
In order that these matrix elements are finite at q 2 — 0, 
the form factors must also satisfy the well known relation 



Ao(0)-^±^^(0) + ^^A 2 (0) = 0. (2) 



2m\ 



2my 



We will work in the static limit of HQET where the 
eigenstates of QCD and HQET Lagrangians are related 
as 



lim 



1 



\H( Ph ))qcd = \H(v)) HQET . (3) 



H 



In this limit it is more convenient to use definitions in 
which the form factors are independent of the heavy me- 
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son mass, namely we propose which in the heavy meson rest frame is the energy of 

the light meson Ey . In such decomposition, all the form 
(V(e v ,p v )\q^Qv\H{v)) = f v e^ uaf: ] e Vv v a p V f), factors (/„, a\, a 2 and a ) scale as constants with the 

(V(e v ,p v )\q^j 5 Q v \H(v)) = -ia 2 (e v ■ v) [p£ - (v ■ p v )v^] heav y meson mass. The relation between the two form 

r *^ , * \ Ml factor decompositions is obtained by corre ctly matching 

~ mi ^ v ~ [V ' €v)V J QCD and HQET at the scale \i ~ m Q HHH: 

-ia (v ■ ey)?A (4) 

where the HQET heavy quark field Q v is independent of 
the heavy quark mass. The form factors /„, a\, a 2 and 
ao are functions of the variable 

m 2 H + my - q 2 , , 

v-pv = » , (5) 

2m H 



[fv{v -Pv) +0{l/m H )\ - 



H run + m v max 



C l0l5 {m Q ) [ao{v . pv) + oil/mH)] = j i ^L_M[2mvA ( q 2 ) + (m H + mv)A 1 {q 2 )-(m H -m v )A 2 (q 2 )] 



V m ff I 9 

+ (m H + £V) Aa ((?2 + my) 2 1 

nifl + my mjy I max 

C7175 i m Q)V™H [a\ (v -pv) +0{l/m H )} = (m H + m v )A 1 (q 2 )\ q 2~ q 2 i ^, 



M{q 



2\ 



17IH + Tfly 



In the following we set the matching constants Cr to 
their tree level values (Cr = 1). At leading order in 
1/roq we thus get 



^l(<7 2 )lg 2 «9^ 
A 0(q 2 )\q 2 ~ql a , 



K fv(v- PV ), 



-ai(v-p v ), 



H 



K a 2 (v-p v ), 



2my 



K a (v-p v ), 



(7) 



which exhibit the usual heavy meson mass scaling laws 
for the semileptonic form factors ^(| . This parametriza- 
tion is especially useful when calculating the form factors 
within HM/^T. The individual contributions of different 
terms in the HM^T Lagrangian to various form factors 
can be easily projected out. 



III. PARAMETRIZATION OF THE FORM 
FACTORS 

Next we propose a general parametrization of the 
heavy to light vector form factors, which takes into ac- 
count all the known scaling and resonance properties of 
the form factors. As already evident from Eq. Q, there 
exist the well known HQET scaling laws in the limit of 
zero recoil |4q|. On the other hand in the large energy 
limit q 2 — > 0, one obtains the following expressions for 
the form factors 123 



V(q 2 )\ q ^ 
^i(<Z 2 )b«o 



Mq 2 )L^o 



TUB + fUy 

m H 
2E V 

run + my 
rnn + my 



U(Ev), 



U(Ey)--^(E v ) 



A (q 2 ) 



= 1- 



£ y (£y)«£||(£y),(8) 



2EyiriH J 

where both universal LEET functions £_l and fy scale 

with the heavy meson mass as m H 3 ^ 2 . These scal ing laws 
were subsequently confirmed by means of SCET III23. 



4 



This is important since the LEET description breaks 
down beyond the tree level due to missing soft gluonic de- 
grees of freedom which are however systematically taken 
into account within SCET. 

The starting point is the vector form factor V, which is 
dominated by the pole at t = m 2 H , when considering the 
part of the phase space that is close to the zero recoil. It 
is very easy to see, that the residuum at that pole scales 
as ~ "if/ 2 with the heavy meson mass Q . For the heavy 
to light transitions this situation is expected to be real- 
ized near the zero recoil where also the HQET scaling JZJ) 
applies. However, since the kinematically accessible re- 
gion q 2 S (0, <?max] i s large, the pole dominance can be 
used only on a small fraction of the phase space, i.e. for 
|q| ss 0. Even in this region the situation for H — > V 
form factors is more complex than in the case of H — > P 
transitions, where g^ ax is indeed very close to the vector 
pole due to low mass of the light pseudoscalar mesons. 
Here, due to larger masses of the light vector mesons, 
9max i s pushed away from the resonance pole and the 
form factor may not be completely saturated by it. For 
the sake of clarity and conciseness we, however, in our 
present study neglect such possible discrepancies and as- 
sume complete saturation of the vector form factor in 
this region by the first physical resonance. On the other 
hand, in the region of large recoils, LEET dictates the 
scaling (jSJ). In the full analogy with the discussion made 
in Refs. |3L l47j. the vector form factor consequently re- 
ceives contributions from two poles and can be written 
as 



V(q 2 ) 



1 - a 



H 



(1 — x)(l — ax) 



(9) 



where x = q 2 /m 2 H * ensures, that the form factor is dom- 
inated by the physical H* pole, while a measures the 
contribution of higher states which are parametrized by 
another effective pole at m 2 s = m H »/a. The parameters 

c' H and a scale with the heavy meson mass as c' H ~ m^^ 2 
and a ~ 1 — oq/ itih to ensure the correct form factor scal- 
ing in both small and large recoil regions. 

An interesting and useful feature one gets from the 
large energy limit is the relation between V and A\ j2^| 



[vV)Mi(? 2 )] 



(m H + my) 2 
2EvrriH 



(10) 



which is valid up to terms oc ^/m 2 H [2jj. This relation 
remains valid even when the leading order corrections due 
to soft gluon exchange are taken into account [26l |2^| . 
When combined with our result jnj, it imposes a single 
pole structure on A\. We can thus continue in the same 
line of argument and write 



^i(9 2 )=4£ 



1-6'; 



(ID 



Here £ = m 2 H /(m H + m v ) 2 is the proportionality fac- 
tor between A% and V from IjlUfl . while b' measures 



the contribution of higher states with spin-parity assign- 
ment 1 + which are parametrized by the effective pole at 
m 2 H i* = mjj, /b' . It can be readily checked that also A%, 
when parametrized in this way, satisfies all the scaling 
constraints. 

Next we parametrize the Aq form factor, which is com- 
pletely independent of all the others so far as it is dom- 
inated by the pseudoscalar pole and is proportional to 
a different universal function in LEET. To satisfy both 
HQET and LEET scaling laws we parametrize it as 



A (q 2 ) 



1 



H 



(l-y)(l-a'j/) : 



(12) 



where y = q 2 /m 2 H ensures the physical pole dom- 

— 1/2 

inance at small recoils. Imposing c" H ~ m H and 
a' ~ 1 — a' /mH preserves all scaling laws, while a' again 
parametrizes the contribution of higher pseudoscalar 



states by an effective pole at m H , 



irr 



2 H /a'. The re- 



semblance to V is obvious and due to the same kind of 
analysis 3] although the parameters appearing in the two 
form factors are completely unrelated. 

Finally for the A2 form factor, due to the pole behavior 
of the A\ form factor on one hand and different HQET 
scaling at <?^ ax (JJJ on the other hand, we have to go 
beyond a simple pole formulation. Thus we impose 



A 2 (q 2 



(m H + m v )£,c' H (l — a) + 2m v c^{\ — a') 



m v )(l-b'x)(l-b"x) 

(13) 

which again satisfies all constraints. Due to the rela- 
tion J2J we only gain one new parameter in this formula- 
tion, b" . This however causes the contribution of the 1 + 
resonances to be shared between the two effective poles 
in this form factor. 

At the end we have parametrized the four H — * V 
vector form factors in terms of the six parameters c' H , a, 
b', a', c'Ij and b". 

It is convenient to introduce hclicity amplitudes for the 
decays H — > Viv as in for example [48| : 



H±(y) 

Mv) 



+ (m H 



■ m v )Ai{m H y) =p ; V{m H y) 



m 2 H {\ - y) 



ran + mv 
2m H mvy/y 

2m H \p v (y)\ 
m v (m H + m v )y/y 



ma + my 
- m^Aximjjy) 



A 2 (m 2 H y) 



(14) 



where y = q 2 /m 2 H and the three-momentum of the light 
vector meson is given by: 



[m 2 H {l-y) 



Am 2 H 



ly. 



(15) 



As shown in Ref. |2j 



these helicity amplitudes can be 
q 2 

relations (JSJ), valid in the large energy limit, one can write 



related to individual form factors near q 2 = 0. Using 



H-(v)\v* 
H+(y)\y* 



2(m H 
0, 



m v )A l (m%y), 



(16) 
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but also, by using relation J3J 



Thus in this region we can probe directly for the param- 
eters c' H (l — a) and d' H (\ — a'). 

On the other hand in the region of small recoil {\pv\ — 
or y ss j/ max ) the helicity amplitudes are saturated by 
the Ai form factor 



H ±(y)\y~y„ 
H o(y)\y~y n 



(m H + m v )A 1 (m 2 H y), 

-2{m H + m v )—A 1 (mj 1 y). (18) 



Consequently we can also directly probe for the value 
of the b' parameter determining the position of the first 
effective axial resonance pole by taking a ratio of H— 
helicity amplitude values at small and large recoils 



H-(y)\ 



y~o 



H -(y)\y~y B 



2 [1 - b'(m H ~ m v ) 2 /m 2 H ] . (19) 



There exists another field multiplet in HMxT contain- 
ing positive parity heavy meson states, — 1/2(1 + 



(17) ma 



-1/37^(7^-^))], where T[ 



is the tensor field with spin-parity assignment 2 + , while 
T2 is another 1 + axial vector meson field. However, as 
pointed out in Ref. , the matrix element of the HM^T 
bosonized currents containing these fields between a sin- 
gle heavy meson state \H(j>h)) and the vacuum vanishes 
at leading order in l/mu due to heavy quark spin sym- 
metry. Consequently, such fields do not contribute at 
leading order to H — > V semileptonic decays. 

Finally we also want to include the radially excited 
states into our discussion and therefore introduce an- 
other odd parity heavy meson multiplet field H' = 
1/2(1+^) [-P M * 7 M — -P'7s] containing the radial excitations 
of ground state pseudoscalar and vector mesons. Such 
excited states were predicted in I The strong inter- 
actions between these fields, ground state heavy meson 
fields H and light vector fields can again be described by 
the lowest order interaction Lagrangian analogous to 122|) 



IV. THE MODEL 



A. Strong interactions 



-iC(H b Vf,p^ a H a ) 
-ip(H b a^ F^{p) ba H a ) 



h.c. 



(23) 



At leading order in chiral and 1 /tuh expansion, strong 
interactions between lowest lying pseudoscalar and vec- 
tor heavy meson fields, and light vector meson fields are 
described by the interaction Lagrangian [TtI ITgj 



-if3{H b v^ ba H a )+i\{H b a^F^(p) ba H a ), (20) 



where the first term is even under parity transforma- 
tion, and the second term is parity odd. H = 1/2(1 + 
j^)[P*7 M — P75] is the matrix representation of the heavy 
meson fields, where P* and P are creation operators for 
heavy vector and pseudoscalar mesons respectively. Light 
vector meson fields are described by p^ — i^=p^ l: where 
p^ is the light vector meson field matrix 



U^ + p1) 



V2 



Pn 



Pi 

K* 



73' 



pi 










Pi) 




(21) 


K*° 









The gauge field tensor is defined as F Al!/ (p) = d^p v — 
dvpfi + [Pn-iPv]- Furthermore, (. . .) indicate a trace over 
spinor matrices and summation over light quark flavor 
indexes. 

In order to incorporate positive parity heavy meson 
states into the model, we introduce the scalar-axial field 
multiplet G = 1/2(1 + ^)[5^7 M 75 - S] representing axial 
(S 1 *) and scalar (S) mesons and incorporate it into the 
interaction Lagrangian by adding additional leading or- 
der interaction terms between heavy even and odd parity 
fields and light vector fields: 



mi -= -K{SbV^G a ) + h.c. 
+ip(H b a fw F^(p) ba G a ) 



h.c. 



(22) 



B. Weak interactions 

For the semileptonic decays the weak Lagrangian can 
be given by the effective current-current Fermi interac- 
tion 



£ eff = -^[^(l- 7 >^] 



(24) 



where Gf is the Fermi constant and J is the effec- 
tive hadronic current. In heavy to light meson decays 
it can be written as J = K a J a , where constants K a 
parametrize the SU (3) flavor mixing, while the lead- 
ing order weak current J a in chiral and l/mji expan- 
sion between heavy ground state pseudoscalar and vector 
mesons and light vector mesons can be written as 113 



+a 1 { 1 h H b pl)+a 2 {Y L l b H b v a pt a ). (25) 

For our leading order calculation, we will also need the 
weak current operator for the scalar and axial heavy 
mesons' transition to the vacuum 



X M = ^'<7 M (l-7 5 )G a ), 



(26) 



and the same for radially excited pseudoscalar and vector 
fields 



(27) 
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V. FORM FACTOR CALCULATION 
A. HMxT calculation at zero recoil 

In HMxPT the derived Feynman rules are valid near 
zero recoil {\pv\ — 0). For the heavy meson propagators 
we use iS a b/2(v-k-A) and -iS a b(g fJ ,u~v, 1 v v )/2(vk- A) 
for the pseudoscalar(scalar) and vector (axial)mesons re- 
spectively, where k^ = q^ L — to#v m . A = is the mass 
splitting between the heavy resonance meson R and the 
ground state strangeless heavy pseudoscalar meson. It 
comes from leading order I /run (spin symmetry break- 
ing), chiral and ST/ (3) breaking corrections, when all 
heavy meson fields are normalized to physical masses of 
ground state strangeless heavy pseudoscalar mesons. For 
the hadronic current matrix element we thus get 

{V{p v )\J"\H{v)) = -iV2g v (ai4 - a 2 v • e y /) 



Xe^v vPVa e vp 



— v2gva \2gva- 



v ■ p v + A H » 



v ■ pv + A H > 



,g v f3v ■ e v v^ . g v ~ C v ' ^v^ 
-i— =oi ; i—=a 



V2 v ■ p v + A Hp \f2 v ■ pv + A H > p 
9v , £y (C - 2 V V -Pv) + {^Pv -(v^V-ey 



-i—=a 
V2 



v ■ p v + A Ha 



(28) 



From this we extract the form factors V(q 2 ), A 1 (q 2 ), 
A 2 (q 2 ) and A (q 2 ): 



v (<l 2 )\q 2 ~ql ax = -^am H Vrn H 

gv 



A 



v ■ p v + A H * 

--=am H ^m H ^— 

V2 v ■ p v + A H >, 

gv , Vm H C ~ 2/w p v 
y/2 m H + m v v -p v + A Ha 



-V2g v ai 



H 



gv , m H + m v n 
-a 



■y/2 V™h v ■ Pv + A H A 



Mq 2 )U 



2V2 m v 



+a 



ft 



c 



v ■ p v + A Hp v ■ p v + A H > 



(29) 



B. Extrapolation to higher recoils 



In order to extrapolate our HMxT calculation results 
at q 2 sa g^ ax to higher recoils we employ the general 
analysis from Sec. III. We model the form factors' q 2 be- 
havior using the formulas ©, (|TT|> . (|T2"|l and (|T3|) with 
model matching conditions at q 2 niix - In order to reduce 
the number of free parameters in this extrapolation we 
employ the same strategy as in our previous work 



We use the information on the contributions of differ- 
ent resonances to the form factors as suggested by our 
model. For the vector form factor V we thus propose 
a = mjj, I m 2 H , t which saturates the effective second pole 
by the first vector radial excitation H'* . Similarly we set 
b' = m 2 H , l'm 2 HA and a 1 = m 2 H /m 2 H , saturating the poles 
of the Aq and A\ form factors and the first pole of the 
A 2 form factor with the H ' pseudoscalar radial excitation 
and the Ha orbital axial excitation respectively. Since 
our model does not contain a second resonance contribu- 
tion to the A 2 form factor, we impose b" — 0, effectively 
sending the second pole mass of this form factor to infin- 
ity. 

At the end we have fixed all the pole parameters ap- 
pearing in the general form factor parametrization for- 
mulas of Sec. Ill using physical information and model 
predictions on the resonances contributing to the various 
form factors. The remaining parameters (ch and c' H ) are 
on the other hand related to the parameters of HMxT via 
the model matching conditions at zero recoil. 



VI. NUMERICAL RESULTS 

We now leave the general discussion of H — > V transi- 
tions and restrict our present study to D decays, although 
our calculations can readily be applied to semileptonic 
decays of B mesons once more experimental information 
becomes available on excited B meson resonances. In 
our numerical analysis we use available experimental in- 
formation and theoretical predictions on charm meson 
resonances. Particularly for the D s axial resonance we 
use the mass m DsJ ( 2i Q ) = 2.459 GeV, while for the first 
orbital excitation of the D* meson, we use the mass of 
TO Di(2420) = 2.422 GeV For the radially excited 

vector resonance we then have the Selex Z)j~j (2632) state 



with mass m 



Dfj(2632) 



2.632 GeV [4fJ. It is important 



to note, however, that so-far the Selex discovery has not 
been confirmed by any other searches [5l|. In D decays, 
the situation is similarly ambiguous. Althou gh t he vec- 
tor D * resonance was discovered by Delphi [53] with a 
mass of m D '* — 2.637 GeV and spin-parity 1~, its exis- 
tence was not confirmed by other searches IH^ 15^. On 
the other hand recent theoretical studies |49tl55l | indicate 
that both radially excited vector states of D as well as 
D s should have slightly larg er masses of m D >, ~ 2.7 GeV 
and m D u ~ 2.8 GeV j4jj. We use these theoretically 
predicted values in our analysis as well as for the radi- 
ally excited pseudoscalar states, for which currently no 
experimental indications exist. 

In our calculations we use for the heavy meson weak 
current coupling a — fH^m*H jESH^i which we derive 
from the lattice QCD value of fo = 0.225 GeV |5al and 
experimental D meson mass mo = 1-87 GeV [5(1 vield- 
mg a — 0.31 GeV 3/2 . The A coupling was usually [HHi 
determined from the value of V(0). However, this deriva- 
tion employed a single pole ansatz for the shape of V(q 2 ). 
One can instead use data on D* — > Dj radiative decays. 
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TABLE I: The pole mesons and the flavor mixing constants 
Khv f° r the D — > V semileptonic decays. 



H V H* Hp H A Khv 



D° 


K*~ 


d'+, d';+ 


Df, D'+ 


D s j(2463)+ 


Vcs 


D+ 


K*° 


d*+, d';+ 


Df, D'+ 


£> sJ (2463)+ 


Vos 


Df 


<f> 


D* s + , D'*+ 


Df, D'+ 


£> sJ (2463)+ 


v ca 


D° 


P~ 


D*+, D*+ 


D+, D'+ 


£>i(2420) 


v cd 


D+ 


P° 


D*+, D*+ 


D+, D'+ 


£>i(2420) 




D+ 




D*+, D* + 


D+, D + 


£>i(2420) 


73 ^ 


Df 




D*+, D'* + 


D+, D + 


£>i(2420) 


V cd 



Following discussion in Refs. [601 l6l| . using the most re- 
cent data on D* radiative and strong decays |50| . and ac- 
counting for the SU(3) flavor symmetry breaking effects, 
we calculate A = —0.526 Gev -1 . The coupling j3 ~ 0.9 
has been estimated in Ref. [§2 relying on the assumption 
that the electromagnetic interactions of the light quark 
within heavy meson are dominated by the exchange of 
p°, lo, <fi vector mesons. 

We fix the other free parameters 
(ai, ait &, Ci A*j Ci A) appearing in the HMxT 
Lagrangian and weak currents by comparing our 
model predictions to known experimental values of 
branching ratios B(D° -> K*~£ + u), B(D S -> (f>£+p), 
B(D+ p°£+u), B(D+ -> K°*l+v), as well as par- 
tial decay width ratios T L /T T (D+ -> K°*i + v) and 
r+/r_(D+ -> K°*e+is) 0- In order to compare the 
results of our approach with experimental values, we 
calculate the decay rates for polarized final light vector 
mesons. Using helicity amplitudes H+ t - t o defined in 
Sec. Ill and by neglecting the lepton masses we get fl7| : 

Fa = GW ^" V? ydy\H a (y)\ 2 \p v (y% 00) 
where a = + , — , and 

V m H J 

The constants Khv parametrize the flavor mixing rele- 
vant to a particular transition, and are given in Tabled 
together with the pole mesons. The transverse, longitu- 
dinal and total decay rates arc then given trivially by 

r T = r+ + r_, 
r L = r„, 

r = r T + r L . (32) 

Consequently, the A form factor does not contribute to 
any decay rate in this approximation and we can not fix 
the parameters «2 and £ solely from comparison with 
experiment. Although A actually does contribute indi- 
rectly through the relation at q 2 = as manifested 



by Eq. 1|17|) , this constraint is not automatically satisfied 
by our model. On the other hand, we can still enforce 
it "by hand" after the extrapolation to q 2 = to ob- 
tain some information on these parameters. Due to the 
specific combinations in which the parameters appear in 
Eqs. (|29(l we are further restrained to determining only 
the products a/2, a'C, and a'/j, using this kind of analysis. 
Lastly, since the only relevant contribution of a\ is to the 
A\ form factor, we cannot disentangle it from the influ- 
ence of a 1 X- Yet again we can impose the large energy 
limit relation l|10(l to extract both values independently. 

We calculate the result for a/2, a'C, a'/i and a>i by a 
weighted average of values obtained from all the mea- 
sured decay rates and their ratios taking into account for 
the experimental uncertainties. Furthermore, the val- 
ues of ot\ and a'C are extracted separately by minimiz- 
ing the fit function {V{Q)£ - A^O)) 2 / (V (0)£ + Ai(0)) 2 . 
Both minimizations are performed in parallel and the 
global minimum is sought on the hypercube of dimen- 
sions [—1, l] 4 in the hyperspace of the fitted parameters. 
At the end we obtain the following values of parameters: 

dtp, = 0.090 GeV 1/2 

a'C = 0.038 GeV 3/2 

a'n = -0.066 GeV 1/2 

ai = -0.128 GeV 1/2 (33) 

These values qualitatively agree with the analysis done in 
Ref. [l8l | using a combination of quark model predictions 
and single pole experimental fits for all the form factors. 

We next use these values in relation (J2J to extract infor- 
mation on the the parameters a2 and C- From Eqs. I|29() 
it is easy to see that the solutions lie on a straight line in 
the a2 x aC plane. We draw these for the various decay 
channels used in our analysis in Fig. ^ We can see that 
all the decay channels considered fit approximately the 
same solution in the plane. Consequently, we can use 
any point on the approximate solution line to obtain the 
same prediction for the q 2 dependence of the Aq form 
factor. 

It is important to note at this point that due to a high 
degree of interplay of the various HMxT parameters in 
the model predictions used in the fit, the values of the 
new model parameters obtained in such a way are very 
volatile to changes in the other inputs to the fit. Fur- 
thermore these are tree level parameter values and may 
in addition be very sensitive to chiral and l/rnu correc- 
tions. Therefore their stated values should be taken cum 
grano salis. However more importantly, the form factor, 
branching ratio and polarization width ratio predictions 
based on this approach are more robust since they are 
insensitive to particular combinations of parameter val- 
ues used, as long as they fit the experimental data. We 
estimate that chiral and heavy quark symmetry break- 
ing corrections could still modify these predictions by as 
much as 30%. 

We are now ready to draw the q 2 dependence of all the 
form factors for the D° — > K~* , D° — > p~ and D s — > <fi 
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transitions. The results are depicted in Figs. [21121 and^J 

Our predictions for the shapes of the various form factors 
can also be summarized using the general formulas 



V(q 2 ) 

Mi 2 ) 

A 2 (q 2 ) 



V(0) 



(1 — x)(l — ax) ' 

A)(Q) 

(l-y)(l-a'y)' 
l -b'x' 



A 2 (0) 



(1 - b'x)(l - b"x)' 



(34) 



where as before x — q 2 /m 2 H » and y = q 2 /m 2 Hp . These 
expressions are actually simplifications of the form fac- 
tor parametrizations © , (|12|) , (|llfl and (|13fl respectively. 
The parameters V(0), A (0), ^i(O), A 2 (0), a, a', b' and 
b" , which we fix by nearest resonance saturation approx- 
imation and HM^T calculation at q^^, are listed in Ta- 
ble |n] for the various decay channels considered. 

Finally, we calculate branching ratios and partial decay 
width ratios for all relevant D — » Vlvg decays. They are 
listed in Table IIIII together with known experimentally 
measured values where we have marked those used in 
the fit of our model parameters. 



VII. CONCLUSION 

We have investigated the q 2 dependence of the heavy to 
light form factors present in the D — > Vivg decays. First 
we have devised a general parametrization of the H — > V 
form factors and then used it to compute the D — > V 
form factors q 2 dependence in the framework of HM^T. 
Although we restrict our discussion to D decays, the im- 
plications of the general form factor parametrization for 
the semileptonic decays of B mesons are obvious, while 
also our HM%T calculations can easily be applied to these 
decays once more experimental data becomes available 
on excited B meson resonances. Furthermore, the exten- 
sion of our parametrization to tensor current form fac- 
tors is trivial using HQET form factor relations |46l |64| , 
which were found to hold even beyond the leading order 
in heavy quark mass expansion [63, [6y| . 

In our calculations we do not include chiral correc- 
tions and 1/niH corrections which might be important 
in charm decays as we already discussed in . As shown 
in Ref. [67j some of these corrections can in fact be ab- 
sorbed into the leading order parameters while a sizable 
number of additional parameters remains undetermined. 
Since our approach is based on a global fit to existing ex- 
perimental data on D — > V semileptonic decays it does 
not seem possible at present to disentangle the effects of 
these new parameters and fix their values. Eventually, 
once more experimental results become available, it will 
be possible to learn more about some combinations of 
these parameters. We note however that while the val- 
ues of HM^T parameters which are used and obtained in 



the fit of our model predictions may indeed be affected by 
chiral and 1/mn corrections, the predicted shapes of the 
form factors are very robust since they are fixed by the 
underlying general form factor parametrization proposed 
in the text and by the masses of the involved resonances. 
Only the overall size of the individual form factors is de- 
termined by HM^T calculation at (^ax an d could thus 
be sensitive to chiral and 1 / m# corrections as well to the 
input parameters of the experimental fit. 

The presence of charm meson resonances in our La- 
grangian affects the values of the form factors at g^ ax and 
induces saturation of the second poles in the parametriza- 
tions of the V(q 2 ) and Ao(q 2 ) form factors by the next 
radial excitations of ■. and mesons respectively. 

The single pole q 2 behavior of the Ai(q 2 ) form factor is 
explained by the presence of a single 1 + state relevant 
to each decay, while in A 2 (q 2 ) in addition to these states 
one might also account for their next radial excitations. 
However, due to the lack of data on their presence we 
assume their masses being much higher than the first 1 + 
states and we neglect their effects. 

We point out that the single pole parametrization of 
all the form factors used in previous experimental studies 
cannot correctly satisfy known HQET and large energy 
limit constraints. In LEET H3 and SCET 29] stud- 
ies the ratio of V(0) and A\{Q) obtained for B — > plv^ 
within these approaches was compared with experimen- 
tal results for D — > K*£ui which were obtained assuming 
such single pole dependence. 

In the current literature there exist a number of studies 
of the q 2 shape of the D — > V form factors. In the lattice 
simulation of Ref. ^(| the scaling behavior of the form 
factors has been properly included, but the pole/dipole 
fits of the form factors they used provide no physical 
information on the q 2 dependence of the form factors 
beyond the leading poles. However, our results for the 
D — > K*£v£ form factors are in good agreement with 
theirs. Similarly the fits done in the quark model calcu- 
lation of [22| to their particular parametrization do not 
differ a lot from our q 2 behavior of the form factors. On 
the other hand, the authors of have used QCD sum 
rules to derive the q 2 dependence of the form factors. As 
presented on Fig.[3[and in TableHTIour results are in good 
agreement with theirs on values of V(0) and Ai(0) while 
our results for the ^(O) are somewhat lower. In Ref. 
the values of the form factors appearing in D p£vg 
decay have been investigated in the same approach and 
they agree well with results of our model calculation and 
extrapolation. 

We hope that the ongoing experimental studies will 
help to shed more light on the shapes of the D — ► V 
form factors. 
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TABLE II: Our model predictions for the parameter values appearing in the general form factor formulas 13411 for the various 
decay channels considered (6" = for all decay modes as explained in the text). 



Decay 


V{0) 




A /n"\ 




a 


/ 

(I 


I 1 



u — > i\ 


n QQ 


119 


U.DZ 


n qi 


\J.D t 


u.oo 


n 7 a 

U. f 4 


D°^p- 


1.05 


1.32 


0.61 


0.31 


0.55 


0.52 


0.69 


D+ -> K°* 


0.99 


1.12 


0.62 


0.31 


0.57 


0.53 


0.74 


D+ 


1.05 


1.32 


0.61 


0.31 


0.55 


0.52 


0.69 


D+ -^lo 


1.05 


1.32 


0.61 


0.31 


0.55 


0.52 


0.69 


D s ^4> 


1.10 


1.02 


0.61 


0.32 


0.57 


0.53 


0.74 


D s - K°* 


1.16 


1.19 


0.60 


0.33 


0.55 


0.52 


0.69 



TABLE III: The branching ratios and partial decay width ratios for the D — » V semileptonic decays. Our model predictions 
and experimental results as explained in the text. 



Decay 


B (model) [%] 


B (Exp.) [%] 


Tl/Tt (model) 


i±/r T (Exp.) 


r+/r_ (model) 


r+/r- (Exp.) 


D° - K-* 


2.2 


2.15 ±0.35 [5j a 


1.14 




0.22 






0.20 


0.194 ± 0.039 ±0.013 [63] 


1.10 




0.13 




D+ -» K°* 


5.6 


5.73 ±0.35 [50j a 


1.13 


1.13 ±0.08 [50] a 


0.22 


0.22 ±0.06 [50] a 




0.25 


0.25 ±0.08 [50] a 


1.10 




0.13 




D+ -^lo 


0.25 


0.17 ±0.06 ±0.01 [g2] 


1.10 




0.13 




D s -Kj> 


2.4 


2.0 ±0.5 [50] a 


1.08 




0.21 




D s -> K°* 


0.22 




1.03 




0.13 





"Values used in the fit of our model parameters 
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FIG. 2: Our model predictions for the q 2 dependence of the 
form factors V(q 2 ) (solid line), Ao(q 2 ) (dashed line), A\(q 2 ) 
(dotted line) and A2(q 2 ) (dash-dotted line) in D° — * K~* 
transition. 




0.5 



I c 

0.2 0.4 0.6 0.8 1 1.2 

q 2 [GeV] 



FIG. 3: Our model predictions for the q 2 dependence of the 
form factors V(q 2 ) (solid line), Ao(q 2 ) (dashed line), Ai(q 2 ) 
(dotted line) and A2(q 2 ) (dash-dotted line) in D° — > p~ tran- 
sition. 
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FIG. 4: Our model predictions for the q 2 dependence of the 
form factors V(q 2 ) (solid line), Ao(q 2 ) (dashed line), A\(q 2 ) 
(dotted line) and A2(q 2 ) (dash-dotted line) in D s — > <f> tran- 
sition. 



